Study of V(25) decays to XJ/if) 



Oh 
< 



M. Ablikim 1 , J. Z. Bai 1 , Y. Ban 10 , J. G. Bian 1 , X. Cai 1 , J. F. Chang 1 , H. F. Chen 16 , H. S. Chen 1 , H. X. Chen 1 , 
J. C. Chen 1 , Jin Chen 1 , Jun Chen 6 , M. L. Chen 1 , Y. B. Chen 1 , S. P. Chi 2 , Y. P. Chu 1 , X. Z. Cui 1 , H. L. Dai 1 , 
Y. S. Dai 18 , Z. Y. Deng 1 , L. Y. Dong 1 , S. X. Du 1 , Z. Z. Du 1 , J. Fang 1 , S. S. Fang 2 , C. D. Fu 1 , H. Y. Fu 1 , 
C. S. Gao 1 , Y. N. Gao 14 , M. Y. Gong 1 , W. X. Gong 1 , S. D. Gu 1 , Y. N. Guo 1 , Y. Q. Guo 1 , Z. J. Guo 15 , 
F. A. Harris 15 , K. L. He 1 , M. He 11 , X. He 1 , Y. K. Heng 1 , H. M. Hu 1 , T. Hu 1 , G. S. Huang 1 * , L. Huang 6 , 
X. P. Huang 1 , X. B. Ji 1 , Q. Y. Jia 10 , C. H. Jiang 1 , X. S. Jiang 1 , D. P. Jin 1 , S. Jin 1 , Y. Jin 1 , Y. F. Lai 1 , F. Li 1 , 
G. Li 1 , H. B. Li 1 *, H. H. Li 1 , J. Li 1 , J. C. Li 1 , Q. J. Li 1 , R. B. Li 1 , R. Y. Li 1 , S. M. Li 1 , W. G. Li 1 , X. L. Li 7 , 
X. Q. Li 9 , X. S. Li 14 , Y. F. Liang 13 , H. B. Liao 5 , C. X. Liu 1 , F. Liu 5 , Fang Liu 16 , H. M. Liu 1 , J. B. Liu 1 , J. P. Liu 17 , 
R. G. Liu 1 , Z. A. Liu 1 , Z. X. Liu 1 , F. Lu 1 , G. R. Lu 4 , J. G. Lu 1 , C. L. Luo 8 , X. L. Luo 1 , F. C. Ma 7 , J. M. Ma 1 , 
L. L. Ma 11 , Q. M. Ma 1 , X. Y. Ma 1 , Z. P. Mao 1 , X. H. Mo 1 , J. Nie 1 , Z. D. Nie 1 , S. L. Olsen 15 , H. P. Peng 16 , 
■ N. D. Qi 1 , C. D. Qian 12 , H. Qin 8 , J. F. Qiu 1 , Z. Y. Ren 1 , G. Rong 1 , L. Y. Shan 1 , L. Shang 1 , D. L. Shen 1 , 

X. Y. Shen 1 , H. Y. Sheng 1 , F. Shi 1 , X. Shi 10 , H. S. Sun 1 , S. S. Sun 16 , Y. Z. Sun 1 , Z. J. Sun 1 , X. Tang 1 , N. Tao 16 , 

Y. R. Tian 14 , G. L. Tong 1 , G. S. Varner 15 , D. Y. Wang 1 , J. Z. Wang 1 , K. Wang 16 , L. Wang 1 , L. S. Wang 1 , 
M. Wang 1 , P. Wang 1 , P. L. Wang 1 , S. Z. Wang 1 , W. F. Wang 1 , Y. F. Wang 1 , Zhe Wang 1 , Z. Wang 1 , Zheng Wang 1 , 

Z. Y. Wang 1 , C. L. Wei 1 , D. H. Wei 3 , N. Wu\ Y. M. Wu 1 , X. M. Xia 1 , X. X. Xie 1 , B. Xin 7 , G. F. Xu 1 , H. Xu 1 , 
Y. Xu 1 , S. T. Xue 1 , M. L. Yan 16 , F. Yang 9 , H. X. Yang 1 , J. Yang 16 , S. D. Yang 1 , Y. X. Yang 3 , M. Ye 1 , M. H. Ye 2 , 
Y. X. Ye 16 , L. H. Yi 6 , Z. Y. Yi\ C. S. Yu\ G. W. Yu 1 , C. Z. Yuan 1 , J. M. Yuan 1 , Y. Yuan 1 , Q. Yue 1 , 
S. L. Zang 1 , Yu. Zeng 1 ,Y. Zeng 6 , B. X. Zhang 1 , B. Y. Zhang 1 , C. C. Zhang 1 , D. H. Zhang 1 , H. Y. Zhang 1 , 
J. Zhang 1 , J. Y. Zhang 1 , J. W. Zhang 1 , L. S. Zhang 1 , Q. J. Zhang 1 , S. Q. Zhang 1 , X. M. Zhang 1 , X. Y. Zhang 11 , 
Y. J. Zhang 10 , Y. Y. Zhang 1 , Yiyun Zhang 13 , Z. P. Zhang 16 , Z. Q. Zhang 4 , D. X. Zhao 1 , J. B. Zhao 1 , 
J. W. Zhao 1 , M. G. Zhao 9 , P. P. Zhao 1 , W. R. Zhao 1 , X. J. Zhao 1 , Y. B. Zhao 1 , Z. G. Zhao 1 *, H. Q. Zheng 10 , 
J. P. Zheng 1 , L. S. Zheng 1 , Z. P. Zheng 1 , X. C. Zhong 1 , B. Q. Zhou 1 , G. M. Zhou 1 , L. Zhou 1 , N. F. Zhou 1 , 
K. J. Zhu 1 , Q. M. Zhu 1 , Y. C. Zhu 1 , Y. S. Zhu 1 , Yingchun Zhu 1 , Z. A. Zhu 1 , B. A. Zhuang 1 , B. S. Zou 1 . 
^- ; (BES Collaboration) 

1 Institute of High Energy Physics, Beijing 100039, People's Republic of China 
f^S ■ 2 China Center for Advanced Science and Technology(CCAST), Beijing 100080, People's Republic of China 

""^ , 3 Guangxi Normal University, Guilin 541004, People's Republic of China 

' 4 Henan Normal University, Xinxiang 453002, People's Republic of China 

^ ■ 5 Huazhong Normal University, Wuhan 430079, People's Republic of China 

| 6 Hunan University, Changsha 410082, People's Republic of China 

• 7 Liaoning University, Shenyang 110036, People's Republic of China 

8 Nanjing Normal University, Nanjing 210097, People's Republic of China 

9 Nankai University, Tianjin 300071, People's Republic of China 

10 Peking University, Beijing 100871, People's Republic of China 
C$ ' 11 Shandong University, Jinan 250100, People's Republic of China 

12 Shanghai Jiaotong University, Shanghai 200030, People's Republic of China 

13 Sichuan University, Chengdu 610064, People's Republic of China 

14 Tsinghua University, Beijing 100084, People's Republic of China 
15 University of Hawaii, Honolulu, Hawaii 96822 

16 University of Science and Technology of China, Hefei 230026, People's Republic of China 
17 Wuhan University, Wuhan 430072, People's Republic of China 
18 Zhejiang University, Hangzhou 310028, People's Republic of China 

* Visiting professor to University of Michigan, Ann Arbor, MI 48109 USA 
< Current address: Purdue University, West Lafayette, Indiana 47907, USA 
* Current address: University of Wisconsin at Madison, Madison WI 5370, USA. 

(Dated: February 7, 2008) 

Using J/ip — » /i + (i~ decays from a sample of approximately 4 x 10 6 ip(2S) events collected with the 
BESI detector, the branching fractions of tp(2S) — > T]J/ip, J/ip, and anything J/ip normalized 
to that of ip(2S) -» tt+tt" J/ip are measured. The results are B(ip(2S) -> r)J/ip)/B(ip(2S) -> 
Tv+TT-J/ip) = 0.098 ± 0.005 ± 0.010, B(ip(2S) -> tv tt J/ip) /B(ip(2S) -» tt+tt" J/ip) = 0.570 ± 
0.009 ± 0.026, and B(ip(2S) -> anything J /ip) / B{ip{2S) -» tt+tt" J/ip) = 1.867 ± 0.026 ± 0.055. 
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PACS numbers: 13.20.Gd, 13.25.Gv, 14.40.Gx 

I. INTRODUCTION 

Transitions of the type ip(2S) — > XJ/ip comprise a large fraction of the total ip(2S) decay width. They include 
exclusive decays where X = t], 7r°, and tttt, as well as the cascade processes ip(2S) JXcO/i/2i XcO/1/2 — * jJ/i 1 - The 
inclusive branching fraction is measured to be B(tp(2S) — > anything J/V>) = (55.7 ± 2.6)%; the contributions from the 
individual sub-processes are less precisely known [1]. 

These branching fractions are important in understanding the hadronic decay dynamics of vector charmonia [2, 
3], since the inclusive hadronic decay branching fraction is calculated by subtracting them from unity. Present 
branching fractions may indicate a possible excess of the tp(2S) hadronic decay rate relative to the "12% rule" 
prediction from J/ip decays. 

Recently the BES experiment reported new measurements with improved precision for B(tp(2S) — ► n°J/ip), 
B{^{2S) - r,J/V>), and Bty(2S) - 7XcO/i/2)5(XcO/i/2 - using decays ^(2S) - 77-W, ■//V' - (where 

^+^~ = or e+e~) [4]. Previous measurements are few and date back to the 1970's and 80's [5-9]. More precise 

measurements are needed. 

In this paper, we report the results of a different technique for measuring branching fractions for the inclusive 
decay ip(2S) — > anything J/V', and the exclusive processes for the cases where X = 77 and X = 7T7t. We reconstruct 
pairs and determine the number of ?p(2S) — > XJ/ip events in our data sample from the J/i/> — ► peak 
in the invariant mass distribution (see Fig. 1). The exclusive branching fractions are determined from fits 

to the distribution of masses recoiling from the J/ip with Monte-Carlo determined distributions for each individual 
channel. We distinguish the ir + ir~ and ir°ir° contributions from simultaneous fits to event samples with and without 
accompanying charged particles. To avoid a number of systematic errors, the channels of interest are normalized to 
the observed number of n + n~J/ip events; we report ratios of the studied branching fractions to that for B(ip(2S) — > 
7r + 7r~ J ftp). This analysis is based on a sample of approximately 4 x 10 6 ip(2S) events obtained with the Beijing 
Spectrometer detector (BESI) [10] at the Beijing Electron-Positron Collider (BEPC). 

II. BESI DETECTOR AND MONTE CARLO SIMULATIONS 

The Beijing Spectrometer, BESI, is a conventional cylindrical magnetic detector that is coaxial with the BEPC 
colliding e + e~ beams. It is described in detail in Ref. [10]. A four-layer central drift chamber (CDC) surrounding 
the beam pipe provides trigger information. Radially outside of the CDC, a forty-layer main drift chamber (MDC) 
provides tracking and energy-loss (dE/dx) information on charged tracks over 85% of the total solid angle. The 
momentum resolution is <r p /p — 1.7%y/l + p 2 (p in GeV/c), and the dE/dx resolution for hadron tracks for this data 
sample is ~ 9%. An array of 48 scintillation counters surrounding the MDC provides measurements of the time-of- 
flight (TOF) of charged tracks with a resolution of <~ 450 ps for hadrons. Outside of the TOF system is a 12 radiation 
length lead-gas barrel shower counter (BSC), operating in self-quenching streamer mode, that measures the positions 
and energies of electrons and photons over 80% of the total solid angle. The energy resolution is ue/E = 22%/y/E (E 
in GeV). Surrounding the BSC is a solenoid magnet that provides a 0.4 Tesla magnetic field in the central tracking 
region of the detector. Three double layers of proportional chambers instrument the magnet flux return (MUID) and 
are used to identify muons with momentum greater than 0.5 GeV/c. 

Monte Carlo simulations are used to determine efficiencies and the expected recoil mass distributions for the 
various processes involved, including "4>(2S) — > 7Xci/2> ^P{2S) — > 7r7rJ/-0, and i(>(2S) — > r/J/ip with J/ip — > as 
well as the background processes e + e~ — > 7/i + /i _ , e + e~ — > ip(2S), ip(2S) (7)/i + /j~, e + e~ — > 27* — > /i + /i~e + e~ 
(/i + /j~/i + /j~). In each decay, angular distributions are generated according to expectations for that process. The 
agreement between the distributions of cos 9^ for data and Monte Carlo has been checked in separate analyses and 
found to be reasonable. Since the default generator for tp(2S) — > mrJ/tp, produces 5*-wave dipion states, while BES 
has measured a small but non-negligible amount of D-wave in the dipion system [11], Monte Carlo events are weighted 
to give the correct angular and mx distributions. 

III. EVENT SELECTION 

Selected events are required to have more than one and less than six charged tracks. 
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A. Muon Selection 

Events must have two identified muon tracks with zero net charge. The fi tracks must satisfy: 

1. 0.5 < < 2.5 GeV/c. Here p^ is the three-momentum of the candidate muon track in the lab. 

2. p^+ > 1.3 or p^- > 1.3 GcV/c or (p^+ > 2.4 GeV/c. This requirement selects events consistent with 
J/ip decay, while rejecting background. 

3. | cos Opl < 0.60. Here M is the laboratory polar angle of the muon. This requirement ensures that muons are 
contained in the MUID system. 

4. cosOp+p- < —0.85. This is the cosine of the angle between the two leptons in the lab. The leptons from this 
decay are almost back to back. 

5. Both tracks must have N hlt > 1, where N hlt is the number of MUID layers with matched hits and ranges from 
to 3. 

6. |iTOF(A i+ ) — £tof(/0| < 4 ns. Here <tof is the time measured by the TOF counters. This requirement removes 
cosmic ray background. 



B. Selection of J/ip events 

For J /ip — > /i + fi~ candidates, the two tracks must satisfy a one constraint kinematic fit to the J/ip mass. 
Shown in Fig. 1 is the invariant mass distribution of the two muons, m w , for J/ip candidates. A clear peak at the 
J/ip mass is evident above background. The distribution of x 2 values from the one-constraint kinematic fits to the 
J/ip — > hypothesis is shown in Fig. 2. The mass recoiling against the J/ip candidates, mx is determined from 

energy and momentum conservation. 



C. Extra track (tt) selection 

In order to distinguish ip(2S) — ► Tr + Tr~J/ip and ip(2S) — ► 7r°7r° J/ip events, separate mx histograms are made 
for events with no additional charged tracks, those with any number of additional charged tracks, and those with two 
or more additional charged tracks. The first histogram and one of the other histograms arc fitted simultaneously [12]. 
To reduce background and improve the quality of the track momentum measurements, events used for this part of 
the analysis are required to have a kinematic fit x 2 < 7. 

Additional charged tracks can originate from ip(2S) — > tt + tt~ J/ip and ip(2S) — > r)J/ip, r\ — > 7r + 7r~7/7r decays, 
and also from gamma conversions and delta rays. Selection criteria are applied to the additional tracks to enhance 
the selection of low energy pion tracks, such as those coming from the process ip(2S) — > TT + n~J/ip, and reject gamma 
conversions and delta rays. 

1. Tracks must have a good helix fit with: 

2. p,,- < 0.5 GeV/c, where p v is the pion momentum; 

3. | cos^l < 0.8, where 6 V is the polar angle of the n in the laboratory system; 

4- p X y„ > 0.08 GeV/c, where p xy7r is the momentum of the pion transverse to the beam direction (this removes 
tracks that circle in the MDC; 

5. \xt E ^ dx \ < 3.0. xi E ^ dX — ^ dE /dx " )m ' a3 ~ { - dE 'j ldx )^p ; where (dE/dx). meas and (dE/dx) exp are the measured and 
expected dE/dx energy losses for pions, respectively, and a is the experimental dE/dx resolution. 

6. For events with more than one additional tracks, cos 6^ < 0.9, where B^-n is the laboratory angle between them. 
This last requirement reduces contamination from misidentified e+e~ pairs from 7 conversions. 

The mx histograms for events with and without additional charged tracks, selected according to the above 
requirements, are shown in Figs. 3 and 4. 
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FIG. 1: Distribution of dimuon invariant mass, m MM , for events that pass the J/ip — > kinematic fit. Dots with error bars 

are data. Also shown is the fit (solid histogram) to the distribution with signal (long dashed histogram) and background (short 
dashed histogram) shapes. 
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FIG. 2: Distribution of \ (data) for events satisfying a one-constraint kinematic fit to J/tp — > /i + . 
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FIG. 3: Fit of the mx distribution events with no additional charged tracks. Shown are the data (points with error bars), the 
component histograms, and the final fit. For the components, the large, long-dash histogram is tp(2S) —> nnj/tp, the narrow, 
dash-dot histogram is 4>(2S) — ► rjj/ip, the broad, short-dashed histogram is tp(2S) — > 7Xci, the broad, hatched histogram is 
tp(2S) — > 7Xc2, and the lowest cross-hatched histogram is the combined e + e~ — > 7/i + ^i~ and e + e~ — > tp(2S), tp(2S) — ► (7)/i + /i~ 
background. The final fit is the solid histogram. 



IV. BACKGROUND STUDIES 

The backgrounds from e+e~ — ► and e+e~ — » tp(2S),?p(2S) — ► (7)^+^^ in the and distribu- 

tions are measured using data. We begin by considering what backgrounds are expected. 



A. Expected backgrounds 

Possible background processes are e + e~ — > 7/i + /i~, e + e~ — > ip(2S),ip(2S) — ► (7)/i + /i~, and e + e~ — ► 
27* — > ^i + /i~e + e~(^i + /j~/j + /i~). The number of background events expected from the simulation of e + e~ — > 27* — > 
^+/i~e+e~ in the m w distribution is four events, which is negligible. The background in the mx distribution with 
no extra charged tracks is even smaller, and the calculated backgrounds for e + e~ — > fi + fi~ fi + /i~ are determined to 
be smaller still. The numbers of two photon background events are found to be entirely negligible, and this process 
is ignored further. 

The Monte Carlo (MC)-determined level of e+e~ — ► -ffi + [i~ background in the m^ [13] distribution and 
the mx distribution with no additional tracks [14] are 556 ± 20, and 192 ± 10 events, respectively. For e+e~ — > 
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FIG. 4: Fit of the mx distribution for events with any number of additional charged tracks. Shown are the data (points with 
error bars), the component histograms, and the final fit (solid histogram). The dashed histogram is ip(2S) — > -k + -k~ J/tp, and 
the hatched histogram is ip(2S) — ► r/J/ip. There is very little evidence for Xci/2- This distribution is composed predominantly 
of V(25) -f tt+tt-J/V'. 

ip(2S),tp(2S) — > ('y)/j, + n~, the expected background in these distributions are 430 ± 67 and 137 ± 22 events, and the 
total backgrounds are 986 ± 70 and 329 ± 24 events, respectively. The result obtained from fitting the background 
level in the m w distribution (Fig. 1) is 1307 ± 56 events which is larger than the MC-determined background level. 
The background shapes and the numbers of background events estimated for the two processes are similar, so the 
histograms are combined in the fitting of the two distributions. 

B. Inclusive Background 

1. Determination of e + e~ — > 7/i + /i~ background from the cos9j/^, distribution. 

The photon from e + e~ — > is typically emitted along the beam direction, producing a dimuon system 

which is along the beam in the opposite direction. The cosine of the angle of the J/ip in the lab, cosflj/^, shows a 
strong peaking near ±1 for simulated e + e~ — > ■jfi + fj,~ events, as shown in Fig. 5b. Some peaking near ±1 is also 
found in the data, as shown for the case with no x 2 requirement in Fig. 5a. By fitting the cosOj/^ distribution for 
data with the e + e~ — ► -ffi + [i~ Monte Carlo distribution plus a distribution to represent the non-peaked events, a total 
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of 839 ± 98 background events is obtained. This can be compared with the predicted e + e — ► 7/i + /i background 
estimate of 556 ± 20 events in the m aa distribution. The difference is 283 ± f 00. 
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FIG. 5: The distribution of the cosine of the angle of the J ftp in the lab, cos 9j/^, for a.) data and b.) simulated e + e~ — -> 7/i + /i~ 
events. The peaks at |cos#j/^| = 1 indicate some e + e~ — > 7/i + /i~ background in the data. No \ 2 requirement is made for 
these plots but m MM < 3.4 GeV/c 2 . 

An additional source of background, not included in the e + e~ — > 7/i + /i~ simulation, is due to "radiative 
return" to the J/^ peak, e+e~ — > jJ/i/j, J/ip — > This background is similar to the signal in the m^+p- 

distribution and would not be part of the background determined from the fit to the m M + M - distribution. However, 
it would be included in the background determined from the cosOj/^ distribution. The difference between the 
background determined from the cos9j/^ distribution and the predicted e+e~ — > -ffi + [i~ is taken as a measure of the 
radiative return. 



2. Total background from the A0 distribution 

Since the [i + and from e+e~ — > -ffi + [i~, e + e~ — > ij){2S),'^{2S) — > and the "radiative return" 

background processes are coplanar, it is possible to determine their level from the A(/> distributions of the data, where 
is the difference between the angles of the /x + and the ^~(±I80°). Fig. 6 shows the Monte Carlo and data 
distributions for A0. The e + e~ — > 7/i + /i~ and e + e~ — > i/)(2S),il>(2S) — > ( r )))i + ^ distributions show a large peak at 
0°, which is not seen in other simulations. A similar peak is seen in the data. The A</> distribution from data is fitted 
using Monte Carlo distributions, including the combined e + e~ — > 7/i + /i~ and e + e~ — > ip(2S),ip(2S) — > (7)/i + /i~ A<ji 
distribution plus a broad distribution to represent the other processes. The fit determines 1476 ± 150 background 
events in the m w distribution. The result is summarized in Table I, along with estimates obtained by adding the 
"radiative return" determined using the cosOj/^ distribution to the predicted and measured background from the 
to^ distribution. The agreement between the various methods to estimate the background is reasonable. 

C. Exclusive Background 

The predicted background from e + e~ — > 7/i + /i~ and e + e~ — > i/)(2S), ip(2S) — > (7)/i + /i~ in the mx dis- 
tribution is smaller than that in the distribution. However "radiative return" which is not included in the 
Monte Carlo generator is not expected to be reduced by the x 2 requirement that is made in going from the m^^- 
distribution to the rax distribution. 

As was done above, the cos9j/^ distribution (but now with a x 2 requirement) can be used to determine the 
amount of e + e~ — > 7/i + /i~ background in the mx distribution. Using a similar calculation, the amount of background 
is determined to be 413 ± 75 events. The measured background is larger than the predicted background, and the 
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TABLE I: Backgrounds in the m MM and nix distributions from e + e~ — > 7/i + ^i~ and e + e~ — > ip(2S),ip(2S) — > (7)/x + /i~ events 
that pass the selection criteria. Predicted results are from Monte Carlo simulations which do not include radiative return. The 
m MM is determined above by fitting background in the m M(J distribution. Results are given with the radiative return determined 
from the difference between the background determined from the cos Oj/^, distribution and the predicted e + e~ — > 7/i + /i~ 
background added to the predicted and m Mfl results. These may be compared to the determinations obtained fitting the 
distributions, which do contain radiative return. 



Distribution 


Predicted + 


m Mfl + 


Acf> 




Meas. Rad. Ret. 


Meas. Rad. Ret. 


Dist. 




f269 ± 122 


1590 ± 140 


1476 ± 150 


m x 


550 ± 80 




636 ± 161 



difference (221 ± 76) is taken as a measurement of the "radiative return". Note that the two determinations of 
"radiative return" are consistent as expected. 

The Acj) distributions with a x 2 requirement may be used, as above, to determine the total background in the 
mx distribution including "radiative return" . The result is 636 ± 161 events. The amount of background when fitting 
the mx distributions will be constrained to this amount. Table I compares the estimates and final backgrounds 
determined from A<j> distributions for both the m w and mx backgrounds. The agreement between the various 
determinations for each distribution is reasonable, and the total backgrounds are quite small. 
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V. FITTING THE MASS DISTRIBUTIONS 



A. Inclusive channel 



To determine the number of inclusive events, the m M + M - distribution is fit with signal and background shapes. 
The signal shape is obtained from real data using the m^+^- distributions from ip(2S) — > anything J/ip (data) 
events with additional charged tracks. These events are primarily due to it + tt~ J/tp. For the background shape, the 
m p.+tj,- distribution obtained by combining the distributions for Monte Carlo e + e~ — > 7/i + /i~ and e + e~ — > ip(2S), 
ip(2S) — > events is used. Figure 1 shows the fit to the to m + m - distribution. The background distribution 

differs somewhat from the data in the high mass region. Because of this, the fit and background determinations 
everywhere in this analysis are restricted to masses below 3.4 GeV/c 2 . The background in the m„+,,- distribution 
(Fig. 1) below 3.4 GeV/c 2 is determined to be 1307 ± 56 events. The number of events in the signal peak is 44498. 



B. Exclusive decays 



To determine the number of exclusive decays and separate tp(2S) — > ir iT J/ip and ijj(2S) — > tt + tt~J/iP events, 
mx histograms for events with and without additional charged tracks, shown in Figs. 3 and 4, are fit simultaneously. 
There are 20818 and 19846 events in the two distributions. Contributions from the Xco are expected to be very small 
[15] and are not included in the fit. The influence of ip(2S) — > n°J/ip is also small, indeed there is no indication 
of such a component in Fig. 3, and this channel is also not included. The mx distributions for Xci, Xc2, and the 
combined background distribution are broad and rather similar in shape, as can be seen in Fig. 3. Since these are 
difficult to distinguish, the Xc2 to Xci ratio is constrained using calculated efficiencies and the PDG world average 
branching fractions for the two processes. In addition, the amount of background is constrained as discussed above. 

In Fig. 4, it is seen that the Xci> Xc2, and rj contributions are small. When fitting, the ratios of these 
components to those in Fig. 3 are also constrained using Monte Carlo determined ratios. Since the amount of 
background is predicted to be very small in this distribution, it is not included in the fit. The numbers of fitted events 
obtained from the simultaneous fits to the mx distributions of Figs. 3 and 4 and corresponding efficiencies are shown 
in Table II. 

As an indication as to how well the Monte Carlo distributions fit the signal, Figs. 7 and 8 show the data with 
all but one shape subtracted, normalized using the fit results, and compared with the Monte Carlo distribution. 
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FIG. 7: Data with two or more additional charged 
tracks with all Monte Carlo distributions, normal- 
ized by the fit result, subtracted except for -ktt (dots 
with error bars), fit with ip(2S) — ► rnvj/ip Monte 
Carlo (solid histogram) plus a polynomial. 



FIG. 8: Data with no additional charged tracks with 
all Monte Carlo distributions, normalized by the fit 
result, subtracted except for r)J/ip (dots with error 
bars), fit with ip{2S) — > r)J/ip Monte Carlo (solid 
histogram) plus a polynomial. 
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TABLE II: Fit results and efficiencies for inclusive and exclusive decays. The exclusive results are for fitting the distributions 
shown in Figs. 3 and 4 simultaneously. The numbers of events, n, and the efficiencies (eff) for the exclusive cases are the results 
for the histogram in Fig. 3. The zero error on the number of \c2 events is because of the x^lx^i constraint. 



CASE 


n 


Sn 


eff 


Anything (m MM ) 


44498 


232 


0.3488 


ir + n~ and 7r u 7r u 


13952 


186 


0.150 


V 


2121 


117 


0.3597 


Xcl 


2793 


62 


0.3678 


Xc2 


1326 





0.3661 



C. Determination of branching ratios 

Ratios of branching fractions, normalizing to the number of ip(2S) — > tt + tt~ J/ip events, i.e. B ^syl^+J-j/^ ) 
are calculated. The advantage of normalizing in this way is that many of the muon selection systematic errors largely 
cancel, as well as the systematic error due to the \ 2 requirement. 

First the correction factor for 7r°7r° events in the tttt mx distribution with additional charged tracks (/i) due 
to gamma conversions and delta rays is determined: 

0.305ei(7r+7r-) 
~~ 0.305ei(7r+7r-) + 0.182ei(7r°7r°) ' 

where ei(7r + 7r~) and ei(7r°7r°) are the efficiencies for additional charged tracks for 7r+7r~ and tt tt events, respectively, 
and 0.305 and 0.182 are the PDG branching fractions for ip(2S) — > tt + tt~ J/ip and ip(2S) — > 7r°7r° J/ip, respectively [1]. 

N^n+n-) = /iiVi(7T7r), 

where Ni(tttt) is the number of tttt events with additional charged tracks. Since ei(7r°7r°) is small, /i is near 1.0: 
fi = 0.967. 

For X = r] or X = Xc, 

B{j>{2S) -> XJ/VQ = £ ra tioi£i(7r+^-)jVo(A-) 
B{i>(2S) -» tt+tt- J/$) eo(A-)/iJVi(7T7r) ' 

where e ra tioi is ratio of the track efficiencies for detecting one tt track for data and Monte Carlo data, N (X) is the 
number of X events in the histogram with no additional tracks, and e (X) is the efficiency for X's in the same sample. 
For the 7r°7r° case: 

B(i/,(2S) -> A J/jj) _ e ratiol ei (7r+7r-)iV (7r7r) e (^+7T-) 
B{ij){2S) -» tt+tt- J/ip) eo(7r 7r )/i7Vi(7T7r) e (7r 7r°) ' 

where Nq (tttt) is the number of events in the tttt peak in the histogram with no additional tracks. 

For the X — > anything case, the m MAt distribution that is fitted to obtain the number of ip(2S) — > anything J/ip 
events, iV(anything J/ip), has no x 2 requirement. Therefore a correction must be made for the effect of this requirement 
since a x 2 requirement was made on the ip(2S) — > tt + tt~ J/ip distribution: 

B(ip(2S) -> anything J/ip) _ ( e ratio i ei (tt+tt") \ / e x 2 cut iV(anything J/tp)\ no x 2 cut \ 
B(ip(2S) -» tt+tt- J/ip) ~ ^(anything J/ip) J \ N^ir+w-J/ip) J 

The efficiencies in the first fraction on the right of the equal sign are determined from Monte Carlo and have a x 2 
requirement. The term e(anything J/ip) is the efficiency for the combination of processes in ip(2S) — > anything J/ip. 
The term e x 2 cut is the efficiency of the x 2 requirement and is determined using data from the ratio of events with 
additional tracks with and without the x 2 requirement. By requiring an extra charged track, a clean sample of events 
with which to determine this ratio is selected. 

The mx analysis is done using the distributions shown in Figs. 3 and 4. It may also be done using the 
distribution for events with more than one additional charged tracks rather than the one with any number of additional 
charged tracks. Each has its advantages. The first case has a higher efficiency and should have a smaller percentage 
error on the efficiency. The second one should have less background from delta rays and gamma conversions. Our 
final results arc determined from the averages of the results for these two cases and are summarized in Table IV. 
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VI. SYSTEMATIC ERRORS 

Systematic errors are summarized in Table III. Some of them are determined by varying selection require- 
ments. Other contributions are determined by turning off the weighting used for ip{2S) — > ttttJ/i/j events, smearing 
the simulated m x by a — 3 MeV to test the effect of changing the mass resolution, and removing constraints. The 
amount of smearing of mx was determined by comparing the width of a Gaussian fit to the 77 peak for data in Fig. 8 
with a Gaussian fit to the corresponding Monte Carlo distribution. The data width is possibly 2 MeV wider (added 
in quadrature) than the Monte Carlo width. 

As shown in Fig. 5b, the cos9j/^ distribution, where Oj/^ is the angle of the J/ip in the lab, shows some 
peaking near ±1.0, which is indicative of background from the radiative process. To check whether this background 
is being handled correctly, the effect of an additional requirement on this angle (| cos0j/^| < 0.9) is determined and 
included in the systematic errors. The background estimate in the mx distribution has also been increased by 30 % 
to determine the systematic error due to the uncertainty in the amount of background. 

In the study done in Ref. [16], the gamma conversion rate in data versus Monte Carlo data was compared, 
and it was found that the rate was lower in data, (4.0 ± 0.15) %, than in Monte Carlo data, 4.5 %. The effect of this 
difference has been studied using two different methods. In the first, the efficiencies for extra charged tracks in neutral 
decays is changed by 4.0/4.5. In the second approach, a much higher track momentum requirement (p xy7r > 0.15 
GeV/c) is made since gamma conversions are predicted by the Monte Carlo simulation to be more important at lower 
momentum. 

This analysis has been done for two cases. The first uses the distribution with any number of additional 
charged tracks, while the second uses the distribution with more than one additional charged tracks. For our final 
results, the averages of the two sets of values are used, and one half the difference between them is included in the 
systematic errors. Also included is an uncertainty of 2 % for the uncertainty in e rat ioi. 

TABLE III: Systematic error summary. The total error uses 1/2 the difference of the fitting results for the greater than 
zero and the greater than one extra charged track cases in quadrature with all the other errors. Here B(00) = B(tp(2S) — ► 
ttVJ/V), B(+-) = B{1>(2S) - tv+tt-J/^), B( V ) = B(^(2S) - r,J/^), B( XcJ ) = B{^{2S) -> 7Xcj)B( X cJ -» 7-W), and 
B(anything) = B(ip(2S) -> anything J/ip). 



Variation 


B(00) 


BM 
B(+-) 


S(Xcl) 
B(+-) 


B(Xc2) 
B(+-) 


B (anything) 
B{+-) 




(%) 


(%) 


(%) 


(%) 


(%) 


1.3 x mx background 


0.26 


0.67 


5.44 


5.44 




p xy „ > 0.08 -» 0.15 GeV/c 


0.88 


0.34 


0.74 


0.51 


0.22 


lower 7 conversion 


1.11 


1.01 


0.74 


1.04 




| cos 0*1 < 0.8 -> 0.75 


0.17 


0.67 


0.25 


0.0 


0.26 


\XdE/dx\ < 3.0 -^< 5.0 
X 2 < 7 ^< 10 


1.93 


0.67 


0.99 


1.04 


0.73 


0.12 


2.35 


0.25 


0.0 


0.03 


| cos 6» M | < 0.6 -> 0.65 


0.29 


0.0 


0.74 


1.04 


0.19 


ymuid > 1 — >> 


0.47 


3.69 


0.99 


1.04 


0.09 


|frroF(A* + ) - *tof(aOI < 4 ->< 5 


0.29 


0.0 


0.25 


0.0 


0.14 


unweight 


1.70 


5.70 


3.47 


3.62 


0.12 


smear mass (3 MeV) 


0.70 


4.36 


0.74 


0.51 


0.42 


remove r\ constraint 


0.18 


0.34 


0.0 


0.0 


0.17 


remove Xc2 constraint 


1.70 


2.35 


29.0 


46.1 


0.45 


I cosflj/^j < 0.9 


0.29 


3.36 


2.23 


2.59 


0.59 


Cratiol 


2.0 


2.0 


2.0 


2.0 


2.0 


> ^> 1 extra tracks 


4.17 


3.05 


3.17 


2.76 


3.48 


m n+M~ background — ► e + e~ — > 7/i + /i~ 


only 








0.14 


Total 


4.58 


9.82 


29.9 


46.8 


2.92 



Also indicated in Table III are the total systematic errors. They are reasonably small except for ip(2S) — > 7Xci 
and ip{2S) — > 7Xc2- The branching ratios are very different when the Xc2 constraint is removed. 

VII. FINAL RESULTS 

The final branching fraction ratios and branching fractions are shown in Table IV, along with the PDG results, 
including their experimental averages and global fit results. For B(tp(2S) — > 7r°7r° J/ijj)/B(ip(2S) — > n + Tr~J/ip), the 
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PDG does not use the previous experimental results and gives no average value. For the other four branching fraction 
ratios, only one measurement exists for each, and Table IV lists the single measurements quoted by the PDG. Our 
results for ^(anything J ftp) / B(tp(2S) — > tt + it~ J ftp) and B(r)J/tp)/B(tp(2S) — > ir + ir~ J/tp) have smaller errors than 
the previous results. As a check, the sum of the exclusive branching fraction ratios in the top of Table IV add to 1.854 
± 0.094, which is in good agreement with the inclusive branching fraction ratio (1.867 ± 0.026), where the errors are 
the fit errors only. 

To determine the branching fractions, the ratios are multiplied by the PDG2002 value for B(tp(2S) — > 
7r + 7r~ J ftp) = (30.5 ± 1.6) %. The agreement for both the ratios of branching fractions and the calculated branching 
fractions using the PDG result for tp(2S) —> B(ir + ir~ J/tp) with the PDG fit results is good. The \c results are high 
compared to the PDG by about one sigma. 

TABLE IV: Final branching ratios and branching fractions. PDG02-exp results are single measurements or averages of mea- 
surements, while PDG02-fit are results of their global fit to many expeimental measurements. For the value marked by 
*, the PDG gives the reciprocal. The BES results in the second half of the table are calculated using the PDG value of 
B^tv = B{i>(2S) -» TT+n- J/tp) = (30.5 ± 1.6)%. 



Case 


This result 


PDG02-exp 


PDG02-fit 


^(anything J /tp)/ B vll 

B(ir°n°J/iP)/B^ 

B{r)J/i>)/B^ 

B(fXci)B(xci - 7J/VO/B™ 
B(-yXc2)B( Xc 2 -> -yJ/1>)/B m 


1.867 ±0.026 ±0.055 
0.570 ± 0.009 ± 0.026 
0.098 ±0.005 ±0.010 
0.126 ±0.003 ±0.038 
0.060 ± 0.000 ± 0.028 


2.016 ±0.150 [17] 

0.091 ±0.021 [8] 
0.085 ± 0.021 [8] 
0.039 ± 0.012 [8] 


1.828 ±0.036* 

0.60 ±0.06 
0.103 ±0.010 
0.087 ± 0.007 
0.042 ± 0.004 


^(anything J/tp) (%) 
B(7r 7r°J/V) (%) 
B( V J/1>) (%) 

B(-y X ci)B(xci -> 7W) (%) 
B(~f X c2)B(xc2 -> 7^) (%) 


56.9 ± 0.8 ± 3.4 
17.4 ± 0.3 ± 1.2 

3.00 ±0.16 ±0.33 
3.9 ±0.1 ±1.2 

1.84 ±0.01 ±0.86 


55 ±7 

2.9 ±0.5 
2.66 ±0.16 
1.20 ±0.13 


55.7 ±2.6 
18.2 ± 1.2 
3.13 ±0.21 
2.66 ±0.15 
1.27 ±0.08 



VIII. SUMMARY 

In this analysis, the m^+^- distribution for candidate J/tp — ► events is fit to determine B(ip(2S) — > 

XJ/ip,X — » anything). Energy - momentum conservation is used to determine the recoil mass, mi, and the mx 
distribution is fit with Monte Carlo determined distributions to determine the number of exclusive decays. All 
processes are normalized to tp(2S) — ► ir + n~ J/tp to reduce systematic errors. Ratios of the branching fractions of 
^(25) -» r)J/tp, n°n°J/tp, and anything J/tp to that of tp(2S) -> n+n-J/tp are measured to be 0.098 ± 0.005 ± 0.010, 
0.570 ± 0.009 ± 0.026, and 1.867 ± 0.026 ± 0.055, respectively. 
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